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Abstract 

In this paper, second post-Newtonian approximation of Einstein-aether theory is obtained by 
Chandrasekhar's approach. Five parameterized post-Newtonian parameters in first post-Newtonian 
approximation are presented after a time transformation and they are identical with previous works, 
in which 7 = 1, /? = 1 and two preferred- frame parameters remain. Meanwhile, in second post- 
Newtonian approximation, a parameter, which represents third order nonlinearity for gravity, is 
zero the same as in general relativity. For an application for future deep space laser ranging 
missions, we reduce the metric coefficients for light propagation in a case of point masses as a 
simplified model of the solar system. The resulting light deflection angle in second post-Newtonian 
approximation poses another constraint on the Einstein-aether theory. 
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I. INTRODUCTION 



Although Einstein's general relativity (GR) has achieved great success both in experi- 
mental tests and in astrophysical applications during the last few decades, the desire to find 
a gravitation theory consistent with quantum theory together with the ever-increasing pre- 
cision of experiments and astrometric observations has urged many "alternative theories" to 
be proposed. Among them, vector-tensor theories are usually investigated in the research of 
preferred frames and violations of Lorentz invariance, due to the existence of a vector field 
in gravity besides a metric tensor g^^. 

Vector fields without a constraint were considered by Will and Nordtvedt for preferred- 
frame theories of gravity in the 1970s They introduced three parameters ai, 0^2 
and 0:3 to describe preferred-frame effects. But, all previous preferred-frame theories were 
ruled out by gravimeter data. After that, Hellings and Nordtvedt investigated a massless 
vector field in addition to the metric field in the solar system experiments and cosmological 
expansion [sl. In 4], a summary of previous works is given and a general action of a vector- 
tensor theory without a constraint is proposed. However, it is shown that an unconstrained 
condition on the vector norm induces instabilities in those theories [s]. 

On the other hand, vector-tensor theories with a potential leading to violations of Lorentz 



invariance were studied by several authors 



A background dynamical tensor field can 



also break Lorentz symmetry. The simplest cases are scalar fields with a non-zero gradient 
and vector fields. A special case is a unit timelike vector field, which is called "Einstein- 
aathe. theory" or "^theory-, tor =hort ^ This theory ha. been rntensively investrgated by 
several authors in the past few years (see [10| for a review) . The linearized theory of ae-theory 



and the propagation of aether waves are studied in In the aspect of post-Newtonian 
(PN) approximation (a weak-field and slow-motion limit), two Eddington-Robertson-Schiff 
parameters, 7 and /5, are obtained by applying the asymptotic weak field limit of spherically 
symmetric static solutions [Q]. One of the parameters due to the preferred frame, 02, is 
calculated in the low-energy effective theory ISj. Then, the parameterized post-Newtonian 



fP 



*N) formalism is used to attain 10 PPN parameters and other constraints are found 



14| . Also, the radiation damping in ae-theory is investigated in |l5l . Il6| . In astrophysics 
application, although it is shown that there is nonexistence of pure aether stars, regular 
perfect fluid star solutions exist with static aether exteriors 171]. In addition, black holes are 



2 



studied in ae-theory and they are found to be ver y c lose to Schwarzschild solution outside 



its horizon and have a spacelike singularity inside [18|]. Furthermore, numerical simulations 
of gravitational collapse in ae-thoery are performed, in which stationary black holes would 
appear as long as the aether coupling constants are not too large [19]. For the properties of 
non-rotating neutron stars in ae-theory, it is shown that it leads to lower maximum neutron 



star masses, as wel 



as larger surface redshifts at a particular mass, for a given nuclear 



equation of state [20|]. Strong field effects on binary system are also considered in ae-theory. 
There exists a one-parameter family with "small-enough" couplings, which passes all current 
observational tests 

To test vector-tensor theories in the solar system by future high-precision experiments, 
the PN approximation of these theories is needed. Therefore FN corrections of equations of 
motion, equations of light, and other relativistic effects should be derived for experiments. 
Two approaches can achieve this task. One is the FFN formalism in first post-Newtonian 
(IFN) approximation proposed by Will and Nordtvedt jj, In this framework, IFN metric 
is parameterized with 10 FFN parameters, and the differences among different theories of 



JJ] , the FFN approach 



gravity are represented by the values of these parameters. In 
is used to attain the IFN approximation of vector-tensor theories in unconstrained and 
constrained cases. Although there are some efforts to extend this formalism to second 



post-Newtonian (2FN) approximation by Nordtvedt and Benacquista in [2^, l23|, |2J] , which 
introduce a lot of parameters to cover various relativistic theories, the ability of such an 
ap pro ach to describe the physical features at the 2FN level is unclear (see a brief comment 
in 251]). In this work, we focus only on ae-theory, therefore we employ a "theory-dependent" 
approach. This approach, which solves the field equations through iteration, is proposed by 



Chandrasekhar |26l. I27j . who obtained IFN and 2FN approximation of GR. 



In what follows, our conventions and notations generally follow those of 28|]. The sig- 
nature of metric is (—,+,+,+). Greek indices take the values from to 3, Latin indices 
take the values from 1 to 3 and repeated indices mean Einstein's summation. Bold letters 
A = (A^) denote spatial vectors. A dot between two spatial vectors, A ■ B, means the 
Euclidean scalar product. 
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II. ACTION AND FIELD EQUATIONS 



In a general tensor-vector theory of gravity, the Lagrangian scalar density involves a 
metric g^i, and a 4- vector field K^. The action defining the theory reads 



.3 



S = 



IGttG 



-gd x + Sm{i>,gij.i.), (1) 



where c is the ultimate speed of the special theory of relativity, G is an a priori gravitational 
constant, g — detlg^,^) < is the determinant of the metric tensor g^^,^, R is the Ricci scalar, 
ip denotes all the matter fields, K'^ = K^Kx and —05°, where 5° is Kronecker S, is the 
asymptotic value of K^. The Lagrange multipher A constrains the vector field K'^ to be 
— 0^. Here, we respect the Einstein equivalence principle so that the matter fields ■0 do not 
interact with the vector field, i.e. the action of matter is the function of ip and g^^, only. 
Variations of g^i, and Kj^ give the field equations 

+ ©2 + + + + = ^T,., (2) 



and 



where 



"(0) "(1) "(2) "(3) "(4) "(5) — l"^^ 



0l°i) — /o-R/ii/ - -^fogt^uR 

+g,unjo - {fo);,u + foRK^K,, (4) 



A 

fj,l/ J fJ,;A-^ 1/ I ^ i-^'- Al/i-^'- -^y 



©2 = fiK,;xK'' + hKx.,,K 
1 

eg) = 2hK\^K(,,,)-2{hK\^K(,),,^ + g,,{h^^^^^ 

-\f29,.K\K^.^^ + f,K,K^K\KP.^^, (6) 

+ [hK'K(,,^)],x - [hK\,K^)],x - [hK^,K^^'\x, (7) 
eg) = 2UK^K'.^K,,^,K,) + UK^K''K^,xK,,, 
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e(5) 

"(1) 
"(2) 
"(3) 
"(4) 



"(5) 



pi 

2 I J,2\ 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



in which ^g{-) = {■)-fivg^'^, f'^ = df^/d^K"^) (/i = 0,1,2,3,4) and parentheses surrounding 
a group of indices mean symmetrization, for example, K(^^.^^ = {l/2){K^.i, + K^.^). The 
energy-momentum tensor T^^ is 

2c dSm{^,g^^) 



(16) 



and T^y = g^^QupT'^'' ■ Following [25|, l29|, ISOj], we define mass, current and stress density as 





(17) 




(18) 




(19) 



Another way to define a involving the PPN parameters 7 and (3 in IPN is 31 1 

a = T°° + 7T'=^' + Ar00(37 - 2/5 - 1)[/ + C(c-^), 



(20) 



where ?7 is the Newtonian potential. Due to 7 = 1 and /3 = 1 in ae-theory (see a PPN 
parameter summary in Sec. IVB4I) . these two definitions are equivalent to each other in IPN 
approximation. It is worth emphasizing that, in these definitions, the matter is described 
by the energy-momentum tensor without specific equation of state. 

Contracting the field equation ([2]) and substituting it into the field equations ([2]) and Q , 
it leads to 

SttG 



R 



-e{K'){T,, + f,,T) 



^ i=i ' 



1=1 



and 



where 6*^*^ = Q^/ilg'^'^, T = T^^g^^" ^ and the couphng functions are 



i=l 



Jo 



and 



ffiu 



fo- 



According to previous works, we assume the action parameters are 

f - 1 

/2 - --pr 



,2' 



/3 
/4 



C3_ 
^2- 



C4 



(21) 
(22) 

(23) 
(24) 

(25) 

(26) 
(27) 
(28) 
(29) 
(30) 



where (/x = 



2,3,4) are constants. When C4 = and A = 0, the action ([T]) reduces to 
the case in y, 4|; When cq = 0, = 1 and A plays a role of the Lagrange multiplier to 
constrain = — 1, the action ([T]) becomes the ae-theory 

mm\ 

15l . llTl . Il8l . l32l | (see 

Table [ID. 

In this paper, we concentrate on ae-theory only, that is, /o = 1, /i = — Ci, /2 = — C2, 
/s = — C3 and /4 = C4. Corresponding field equations of ae-theory are simplified as 
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R„. = ^ I T,„. - ;(7„,.T 1 + ^Q„,„ V e« - V ew, (31) 



^ ^ i=l 1=1 



and 



(32) 



i=l 



TABLE I: Special cases of the action ([T]) 



Special case 



parameters 



General Relativity 
Einstein-Maxwell theory 
Will & Nordtvedt []J 
Hellings & Nordtvedt \3\ 
Will [4] 

Aether theory [9, 11, 12, U, 15, 17, 18, 32] 



Co = Ci = C2 = C3 = C4 = 
Ci + C3 = 0, Co = C2 = C4 = 
Co = C2 = C3 = C4 = 

ci + C2 + ca = 0, C4 = 

C4 = 0" 

Co = 0, K^ = -1 



"■Co = -uj, ci = 2e - r, C2 = -ij, C3 = t] - 2e 12]. 



and the Lagrange multiplier A from Eqs. (fT5|) and (l32l) can be expressed as 



(33) 



where i^T^ = — 1 is used. These field equations coincide with previous works' 



15 



17 



18 



32]. 
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III. PN EXPANSION OF THE METRIC AND VECTOR FIELD 



In PN approximation, we consider an asymptotically fiat spacetime, whose metric g^i^ to 
second order has the form as 



5'oo 



(2) (4) (6) 

-1 + hoo + hoo + hoo + 0{e'), 



(3) (5) 



9ij = Sij + e hij + e hij + 0{e^ 



where e = 1/c. Furthermore, we simplify the notations with the definitions 331 ]: 



(2) 

N= /ioo, 



_ (6) 
Q = hoo, 



And the expansions of vector field are 



_ (3) _ (2) 

Li = h Qi, Hj^ = h 



_ (5) 
Qi= hoi, 



Qij — h ij. 



_ (2) 

H = hkk, 



(2) (4) (6) 



(34) 
(35) 
(36) 

(37) 
(38) 

(39) 



(3) (5) 



which is a timehke unit vector. Hence, 



(2) 



2 



(2) (2) 



(4) 



— e 



(2) 



(4) 



(2) 



+ 2NL - LkLk + Q + NKq - 2NKq - 2(N^ + L)Ko 



(3) (3) (3) (2) (4) (6) ■ 

+ 2LkKk - KkKk + 2KoKo - 2Ko 



(2) (4) (6) 

Consequently, with the unitary and timehke condition K = —1, Kq, Kq and Kq 



solved at corresponding orders as 



(2) TV 
2"' 



and 



(4) L 



(6) Ar3 LN Q If \ f 



IV. GAUGE CONDITION 



The gauge condition we use for the metric is the harmonic gauge 



99n,u = 0, 



which reads 



and 



-H, 



-Ni — Hikh 



4 1 1 1 
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— -j^HikHik^i — -L^i + Li^t — Qik,k + -^Qkk, 



0{e% 



+ -^LkH^k — LiHik^k — -^L^t — -^Qkk,t + Qk,k^ 

= 0{e'). (47) 
In addition, the divergence of the vector field equation (132!) . which is 

5 



i=l 



a 



can be apphed to simphfy our mathematical deduction whose expression will be given 
in Sec. IV B 21 in detail. 

V. SECOND ORDER POST-NEWTONIAN APPROXIMATION 
A. Newtonian limit 

The leading terms of goo and Kq show the Newtonian limit of ae-theory. With Eqs. (l33l) 
and (12]), Roo to the order 0{e'^) yields 

AN = -SnQa, (49) 

where A = is the Laplace operator for the space coordinates and Newton's constant 
Q is related to the constant G by 

- (50) 



1 — |ci4 

Here, we use a notation like Cijk for q + Cj + Ck, for example C123 = ci + C2 + C3. Obviously, 
is twice of usual Newtonian potential U, which is given by the standard Poisson integral 

U = A-'{-4nga} [ p^d'x'. (51) 

J \x — X \ 

B. First order post-Newtonian approximation 



Following Chandrasekhar's approach 26|, |27|] , we look for the solution of the field equa- 



tions in the form of Taylor expansion with respect to the parameter e. The solutions of the 
metric Qfj,,^ and the vector filed follows. 
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1. Hij 

From the field equations of Rij witli harmonic gauge, we can easily have 

AHij = -8nga5ij. (52) 

As Hij is solved in an isotropic form, it brings a lot of convenience into subsequent works. 
The harmonic gauge and the covariant divergence of the field equation for the vector field, 
Eq. fHHl) . become quite simple in IPN approximation. 



(3) 

2. Li and Ki 

Expanding the field equation of R^i and Ki to 0{e^) with gauge ( 147|) . we obtain 

(3) (3) 

(1 - ci3)ALi + cisAKj - ci23A^,it + (c2 + ci23)Kk,ki = IQirGcTi, (53) 

and 

1 (3) (3) 1 

- -cigALi + ciAKi + C2zKk,ki - -(C23 - Ci)N^it = 0. (54) 

(3) (3) 

To solve ALi and AKi in above equations, we use Eq. (HHj) to eliminate the terms of Kf^^^i. 
With the help of harmonic gauge and the results obtained previously, Eq. (j48l) can be 
written down as 



' 2 2ci23 ; ' . 



O(e^). (55) 



(3) 

Therefore, we can solve L,- and Kj as 



and 



lQc,{2-cu) 2(ci + ciC4) , . 

^"-^ = c^ + 2c,-cr^ ^ + 4 + 2c,-cf^''^ ^'^^ 



(3) 

where Cj^^^ is a constant. All constants with the form of Cy ^ are given in Appendix Rl 



3. L 



As above, we expand the field equation of i?oo to 0{e ) and solve L as 

AL = -iAiV2 + Ci^_^^Ar,,i. (58) 
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4- Summary of PPN parameters 

In IPN approximation, the metric (1MI) - (136I) is equivalent to the PPN metric [l, 4| under 
a trivial gauge transformation. The transformation between our reference system (t, x*) and 
the PPN reference system in the standard PN gauge (tpN, ^^^pn) reads 



where 



Ai 



1 — |ci4 

and X is the superpotential defined by 



7 9 3 1 , 

1 ~ oCl - 7TC2 - T^Ca - 2C4 + (C2 + Ci2 + C1234J 

III /C123 



y a(a;',t)|a;-a;'|dV + 0(e2), 



so that 



Ax = 



After transformation, 5 PPN parameters are 

7 = 1, 
/5 = 1, 
i = 0, 

8 

tti = 

and 



8(CiC4 + C3) 



^2 



(2ci3 - Cm)' I2C3C13 + 2cici4(l - 2ci4) + (c? - cj) (4 - 6C13 + Tcm) 



(59) 
(60) 

(61) 

(62) 
(63) 

(64) 
(65) 
(66) 
(67) 



Cl23(2-Ci4) (2 - Ci4)(c| + 2ci - cf) 

Other 5 conservation law parameters, 03, Ci, C2, Cs si'iid ^4, are all zero due to the theorem 
of Lee, et al. 



34| . These results perfectly match previous works (l2|, 1^ . 



C. Second order post-Newtonian approximation 



Following above procedures, we can obtain 2PN approximation of ae-theory. The 0{e^) 
term in gij is solved as 

/ 1 \ / (3) (3) \ 

^Qij = - 1 1 - -ci4 1 N^iNj + ci3 1 Lijt + Lj^it - Kijt - Kj^it j 



11 



-8(2 - ci4)7r^(Tij + 5i. 



1 



+ ^AiV^ + 8(2 - c^^)7:gakk + C^^N^tt 



(69) 



As we can see, Qij no longer keeps isotropic as Hij. To succeed in the convention and 
convenience in IPN, we will try to transfer Qij into an isotropic form as possible as we can 
in next section for the application in the light propagation model. 

Before solving Qi and Ki, we need Eq. to attain to 0{t ). It gives 



■(3) 



2 

1 /I 



3 Ci2 + C24 \ 



2Ci23 / ' . 



3 1 3 , ^ 

, Ci Co C3 H — C4 — 1 

C123 V2 4 ' 4 ' 4 ^ ' 



NN^t + — ( |ci3 + Jc2 ) NLk^k 
C123 V" ^ 



+ -LfcAT^fc - ( 1 



C4 



(5) 

+ irfc,fc + A 



0{ 



4Ci23 
7\ 



(3) 



2c 



123 



2 VC123 



(70) 



where A^ satisfies a Poisson's equation given in Appendix El After eliminating the terms 

(5) (5) 

related to K^^k in the field equations, we can solve Qi and Ki as 



AQi 



and 



(5) 

A/sT, 



TiQNcJi 



nQaLi 



16ci(2 — Cm) 
c| — cf + 2ci 

^ 16(ciC4 + ci) ^^^g)^ + ,,Ar^,Ar, + C^V .iVAT,, 

C3 — C]^ + ZCi 

2ci 

+ ^5 o . ^ — N^ikLk — 2N^kLk,i H — 2 9 , — N^kLi^k 



{cl + cl + 2ciC4 - 2ci) 
C3 — c^ + 2ci 



NN,u 
C13C14 



- cf + 2ci 
2ci + C13C14 
c| — cf + 2ci 
C13C14 



cf + 2ci 



LkN^ki + Li^tt + 



2ciCi4 



C3 — cf + 2ci 



cf + 2ci 



KkN^ki + 2'^^ \ , Kk^N^k 



2ci3Ci4 C 
Ai,ft + -2 2 



C3 — '^1 2ci 
12ci3(2 - C14) 



C3 — '^i ~l~ 2ci 



cf + 2cr"""" c^-cf + 2ci 

^ _ 2(C3C23 - C1C12) 
^kk,it 2 2 I o "^1*5 

C3 — ^1 ~r 



(71) 



c\ — c\ + 2ci 



, 8ci3(2 - C14) 



cf + 2ci 



4(cf + 2ciC34 + 4 + 2C3C4 - 4C3 



2ci 



C3 — cf + 2ci 



I (3) 

-TcQaKi 



(5) 



(5) 



Cl3 



c^ — cf + 2ci 



-NikLk 



-.NuLki + 



Cl4(Cl3 - 1) 



2" ■ ci-cf + 2ci 



NhLik 



Cf3 + 2C13(C4 + 1) - 2C4 

2(ci - cf + 2ci) 



UN 
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2(ci-c2 + 2ci) 

Cl3(2Ci + C4-C3)-2ci(3) 

ci - c? + 2ci 

^13(^3 — — 2C23) 

"2Ci23(ci-c2 + 2Ci) 



2(ci-c2 + 2ci) 

2C14(C13-I) (g 
c2 - + 2Ci 



i.tt 



- C3 + 2C23 . 
C3 — ^1 ~r -^^1 



Similarly, Q can be solved as 



16(4 + cl + 2ciC4 - 2ci) 
c| — + 2ci 



-KQLkCTk 



Cg — C-^ ~\~ ZC\ ' 



kJ^,k 



+N^kQki,i + N^kiQki — 7^N^iQkk,i 7^ —Lk,iLi^k 



2(1 + C3-C4) 1 

H J-'l,kJ-il,k + 



2 - Ci4 



(3) 



2 - Ci4 

2(ci + 2c3 - C4) 



2 + 



Cl4(Cl2 + C24) 



C123 



+C*^ N^ktKk — [ Lk,iKi^k + Li,kKi,k 



(3) 

NktKk 



(3) 



(3) 



2(C23 - C4) ^ ^ 
-J-'kJ^l,lk + 



2 - Ci4 

2(C1-C4)(3) (3) 



2 - Ci4 
Cl3 



2 - Ci4 



4C3 (3) 



2 - Ci4 



Ki^kKi^k 



2 - Ci4 



T Ci2 + C24\^ , 2(ci2 + C24 + C123) . 

-L H \Uikk,tt H >A,t. 



C123 



2 - Ci4 



(72) 



(73) 



D. Verification of the solutions of metric and vector field 



One way to verify the solutions obtained above is to check the gauge condition. Inserting 
the metric coefficients into the harmonic gauge, we have 

1 



^ ^3 32ci7rG(l + e^iV) 
c| + 2ci — c\ 



and 



<^,t + <7k,k + e ( 2'^-^.* ~ 



1 



(74) 



(75) 



which are equivalent to the equation of motion T.^'^ = 0. 

Another approach as the most reliable way to verify the results is to substitute the metric 
and vector 2PN expansion coefficients into the field equations and check. Our results pass 
this examination. 
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E. The 2PN parameters i 



The 2PN metric of ae-theory shows a picture of second order PPN (2PPN) formahsm. 



22,123 



2J]. Whereas 



Previous works of 2PPN formahsm focus on a many-body Lagrangian 
we have not restricted the matter in our model, we are not going to compare our resuhs with 
theirs in this paper. Despite of applying a "theory-dependent" approach, several works also 
obtain some 2PN parameters. Although there are many parameters regarding the effects of 
preferred frame, here we concentrate only on the parameter that represents the third order 
nonlinearity of Newtonian potential in g^Q only here. Damour and Esposito-Farese calculate 
two 2PN parameters e and C in a multiscalar-tensor theory, and they find that the possible 



2PN deviations from GR, Sg^^ = g^^ — g^^, are given by [25|, |30| 

Aisfo-^ = ^^U' - "-M^oU- + O (2, ^) + O , (76) 

where only one scalar field is involved, U is the Newtonian potential, e measures how much 
the third order nonlinearity there is in the superposition law for gravity and C is no longer an 
independent 2PN parameter with = ({^,(3). Recently, Xie et a/. [35,] find a 2PN parameter 
t in a scalar-tensor theory (STT) with a intermediate range force, in which l in Sg^J"^ = 
9oo'^ — g^i^ has the form as 

A^Jo^ = -^f/Af/Vo(l f) +o(i). (77) 

With the help of AU = -A-nGa and AU^ = WAU'^ - SW^AU, l and e are equivalent and 
they represent the third order nonlinearity in 2PN g^^, which is totally different from the 
terms due to the combinations of IPN terms, 0{pf / , [3 / c^). In the case of ae-theory, t = 0, 
which makes ae-theory no differences with GR in the parameters of 7, (3 and t (see Tab. 
HTj) . In this table, Wq, ^1 and UJ2 are constants, coming from the expansion of the coupling 
function in the STT with an intermediate-range force (see {ssl for details). But even so, 
experiments, especially the deep space laser ranging missions, still can test ae-theory by its 
unique effects that deviate from GR, such as the 2PN light deflection angle. 

VI. A 2PN METRIC FOR LIGHT PROPAGATION IN THE SOLAR SYSTEM 

Future deep space laser ranging missions such as Laser Astrometric Test of Relativity 



(LATOR) [36|, and Astrodynamics Space Test of Relativity (ASTROD) [37|, together with 
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TABLE II: Summary of the parameters. 



Parameter 


What it measures, relative to GR 


Value in GR 


Value in STT [35] 


Value in ae-theory 




How much space curvature 

is produced by unit rest mass? [28] 

How much the second order nonhnearity 

is there in the superposition law 

for gravity (500)? [28] 


1 
1 


<^o+l 
uJo+2 

1 1 

(2LJo+3)(2tJo+4)2 


1 
1 


t 


How much the third order nonlinearity 
is there in the superposition law 
for gravity (500)? 





2{3+2a;o)(aJo+2)^ 






astrometry missions such as Global Astrometric Interferometer for Astrophysics (GAIA) 38 1 
and Space Interferometry Mission (SIM) 39| will be able to test relativistic gravity to an 
unprecedented level of accuracy in the solar system. Those missions will enable us to test 
relativistic gravity to 10~^ — 10~^, and will require 2PN approximation of relevant theories 
of gravity, including metric coefficients, equations of motion and equations of light ray. 

Hence, in this section, we discuss a 2PN metric for light propagation in the solar system. 
Considering a practical model, we just study a situation of N point masses in a global frame 
as the first step. And we impose a constraint on the metric, that is, after ignoring all the 
planets in the solar system, the spatial part of the metric, Qij, should be isotropic after a 
coordinate transformation. 



A. A coordinate transformation 

A coordinate transformation 

= x^' + e^e{x''), (78) 

and 

= x'^-e4^/^(x"), (79) 
where C,'^ ~ ^(1)5 changes the metric to 
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(80) 



When the time component is chosen to be fixed {t = t), we obtain the metric to 2PN order 



= 9^A^n - e%{x-) - e%{x-) + 0{c 
9o^ = go^{xn-e'Ct{^'') + 0{e'), 
900 = ^oo(x") + 0(e^). 



(81) 
(82) 
(83) 



Then, in the right hand side of above three equations need to be replaced with x". We 
must also transform the functional integrals over x'^' that appear in (yf^j, into integrals over 
x'^' . The only place where this changes anything is in goQ = —1 + e^N{x,t) + O(e^), where 



N{x,t) = 2g j 



^^d'x'. 



\x — x'\ 

Like Eq. (!80l) . the transferred energy-momentum tensor T^'^ is 

_ t^x^ dx" 
dxP dx^ 

= T^'' + e^T^P^p + t^T^'Pi^p + O(e^). 
According to the definition of a f|T7j) . it can be obtained that 



(84) 



(85) 



- = rpOO^rpSS 

rpOO _j_ rpSS _j_ 2,^^T^^ 

= a + 0(e^). 



(86) 



Furthermore, the difference between the volume element d^x' and d^x' is a Jacobian deter- 
minant that 

dx' 



d'x' 

where $, = ^{x). We also have 

1 



d^x' 



dx' 

dV[l-e"V'-|' + 0(e«)]. 



(87) 



\x — X 



\x-x' - - I' 
1 



\X — X' 



+ e 



,(^-x')-(^-O^0(^8)_ 



\X — X 



'13 



(88) 
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Thus, we put these relations together and have 

N(x,t) = N(x,t)-2e^g I ^dV 

J \x ~ x'\ 



J \x — x'[^ 

Finally, we obtain the metric after transformation 

= 9^A^n - e%{xn - e%{xn + 0{e'), (90) 
9o^ = 9o^{xn-e'^U^n + 0{e'), (91) 
^00 = ^oo(x")-2e66?A-H-47raV-|} 

+2e'g [ ^^i^^l^^Sx' + 0{e'). (92) 
J \x — x']-^ 

(3) _ (3) (3) 

So if = C13 A {Li^t — Ki^t), the term of C13 A {Li^t + Lj^ — Kijt — Kj^) can be eliminated 
and the IPN parts of metric keep unchanged. With the help of this transformation, the 
metric for light can be changed into a simpler form. 



B. A 2PN light propagation metric of N point masses 

In light propagation, we can cut off the full 2PN metric to 

goo = -l + e^N + t^L + 0{e'') (93) 
go^ = e^L, + 0{e'), (94) 
g,j = + e^Hi, + e^Q,, + 0{e'). (95) 

em of non-spinning point masses for 



If we consider the solar system A^-body prob 
simplicity, we follow the notation adopted by [4^, 



4l| and use the matter stress-energy 



tensor 

c'T>^^{x,t) = Y,f^a{tKitKm^ - ya{t)), (96) 

a 

where 6 denotes the three-dimensional Dirac distribution, the trajectory of the ath mass 
is represented by Vait), the coordinate velocity of the ath body are Va = dya(t)/dt and 
= (c, Va) and /Xq denotes an effective time-dependent mass of the ath body defined by 



'-y^99pxVa 
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'-'a ^ a 



where means evaluation at the ath body and being the constant Schwarzschild mass. 
Another useful notation is 

Ut)=^^a{m + t\l), (98) 

where v"^ = v"^. Both /i^ and fla reduce to the Schwarzschild mass at Newtonian order: 
fJ'a = + 0{e^) and /ta = rria + 0{e^). Then the mass, current and stress densities ( |T71fT9|l 
for the point masses read 



a 

a 
a 

Therefore, we can work out and Hij quickly. 



A^ = 2A-\-A7Tga} 



82 \^ Qfrib 



by the relation of /ia that 



/ia = rna< 1 + e 



rrial l + t 



N--h\ +-vl 
2 L 2 " 



2 j-" Tab 

bfa 



where = |a? - ya\ and Tab = \ya - Vb] 
In the solution of L,-, due to 



A-'i-ATTGa,} = E / 



gmavl6{z - y 



\x — z\ 



GrUaVl 



(99) 

(100) 
(101) 



(102) 
(103) 



(104) 



(105) 



and 



= ^{^SYl j ma5{z-ya)\x-zlA^z + 0{^] 
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Q-rria 



where 



it can be solved as 



X^^g J a{x',t)\x-x'\d^x' + 0{e'), 



where Ua ^ {x - Va) / Va- 
In qqq, with the help of 



\X - ya\,t = -ria ■ Va 



and 



we can obtain 



X - ya\,tt = K 2^ • nab 



= a|^;^ J mj{z-ya)\x-z\d'z + 0{e')^ 
= AlgJ2ma\x-ya\,tt + 0{e^)\ 



grria 



vl - {riaVaf 



which leads to 



i^--2E^+c*,..ET^ 

' n I n. 



-2 + ^Nu 2^ 2^ —r2 ■ riab + 0{e), 



ran 



a by^a 

where Hab = {ya - yb)/rab- 

The quadratic part of potentials in Q^j can be rewritten as 

-.-.^4Eex(i)/i)^^.4EE^w.(i)/i 



(106) 



(107) 



h = ^ ^ + f^tf;^ E ^(n„ . .„)< + Oie% (108) 

c| - c| + 2ci ^ r„ c§ - cf + 2ci ^ ra 



(109) 
(110) 



^J2J2^^^na-nab + 0{e')i (111) 



(112) 
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\ E ^'-'(^5 + ^^A) (;^) + 4 E E Q'rn^rn,d^,d,, (J-) , (113) 



where dai denotes the partial derivative with respect to y^. The integral of the self-terms 
can be readily deduced from A(lnra) = 1/r^ on the other hand, the interaction terms 



are obtained by 



AlnS 



ab 



ran 



;ii4) 



where Sab = ra + rf, + Tab |42|]. Consequently the first term in AQjj can be solved as 

^ ^ a b^a 



a a 



+4 E E ^^"^'^"^'^ 

a b^a 



"Tab Sab 



C2 
'-'ab 



(115) 



where nj|, = n'^ifn^ and two relations that 



6,, - 2nyi 

L/,, 111 la o ) 



;ii6) 



and 



are used. 



OaiObj m bab = H 



rabS, 



ab 



C2 



(117) 



(3) 



In solving Qij, we can use the transformation fl90l) - fl92l) and set = C13A ^(Lj ^ — Ki^t) 

(3) (3) 

to eliminate the term C13A {Li jt + Lj^n — Kijt — Kjn). After that, with the relations that 

A-'{-4nga.,} = E / ^^^Siz - ya)d'z 



E 



118) 



and 



(119) 



transferred qI^^ can be worked out as 



Q 



(1) 



+4 1- 



E 



1 \\^r2 
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4M - ici4 ) ^^^^mamft 



a b^a 



^ab ^ (^4 - K) «fe + nl) 



Grria 



fab Sab 
,2 



C2 



4 1- 



ran 



a b^a 



Collecting all these results together, we have the metric for 2PN light propagation as 
2gma 



E 



ra 



+ E E g^rriamb 

a by^a 



2 



2E 



TaTb TaTab 



2 ■ Tlab 



' ab 



i/Oi 



,3 J ^"-IS 



2c\^ - 2ciCi4 + 8ci >^ Qruavl , 2(03 + C1C4) ^ Qrria , 



C3 — '-'1 ~l~ 



2^ma 



E 



+ 



C3 — '-'1 ~l~ 



;i2i) 



(122) 



Qnia 



a b^a 







^rriamb 


2 2 


TaTb TaTab 



\ / ^ a 

Ua ■ riab , 2(2 - Ci4 



ab 



+^'{ +4(1 - E^.x+ (1 - 



2_2 ^a^a 



-4 1 



-5=")EE5^ 



marrib 



a b^a 



^ab , Kfe-<)K6 + ^6) 



fabSa 



+ 



C2 
•^ab 



+0(e^). 



(123) 



When ci = C2 = C3 = C4 = 0, ae-theory goes back to GR, and the above metric reduces to 
the metric given in GR by *4u] . 

If we only consider the case that light just passes the limb of the Sun (Mq), which will 
provide the strongest light deflection effect in the solar system, we can simplify the 2PN 
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metric further, by neglecting the contributions from planets. Hence, the metric for such a 
practical 2PN light deflection experiments can be simplified as 



(2) 
^00 

(2) 
(2) 

9ij 



Rc. 



26" 



& 



-1 + e- 
0, 



R^ 



-Oi-i 



■& 



+6" 



1 + -cu 



Rq 



1 ^ -2 nj2 ^0^0 



(124) 
(125) 

(126) 



which is anisotropic in Qij, causing unconvenience in the calculation of 2PN light deflection 
angle. Only one 2PN parameter Ci4, which deviates from GR, remains in this special case. 



C. Isotropic coordinates for the dominated body 



The transformation between isotropic and harmonic coordinate system in GR involves 
only the radial coordinate, and is given by 



th + ^ ym = ri\ 1 + e - — 

2rj 



whose 2PN approximation is 



n = rH- e 



4r 



H 



(127) 



;i28) 



where r/ and th represent the radial coordinate in the isotropic and harmonic system re- 
spectively. Hence, we obtain the transformation in ae-theory for spatial components as 



1 \ g^m^ 



a;'-eNl--ci4 



-n 



(129) 



2 J 4r 

and, using the transformation fl90l) - fl92l) again, we get the metric g^iJ for 2PN light propa- 
gation as 

2gMQ 2g^M^ 



c^r 



+ 



3,1^ \g^Mi 



\^2gM^ 



(dV^ + r^dn^) 



(130) 



where df2^ = d^^ + sin^ ^d0^. In the spacetime with metric (1130p . the light deflection angle 
up to 2PN approximation is 



'15 + Ci4)7r 



2 ]\/r2 





g'Mi 



:i3ii 
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where d represents the coordinate radius at the point of closest approach of the ray and 
Ci4 is the only non-GR parameter representing the deviation from GR, which poses another 
constraint on action parameters of ae-theory. 

Many previous works have been done to constrain the parameters Ci^2,3,4 by theoretical 
and experimental analyses, including the rate of primordial nucleosynthesis 43|], the rate 
of Cerenkov radiation [51], the stability and positive energy of the linearized wave modes 



HQ, 



the experimental bounds on the PPN parameters 



13 



14l |. which shows a two- 
parameter family of ae-theory can satisfy the requirements of observations (see for a 
summary). The constraints coming from the strong field effects and the rate of gravitational 



radiation damping in binary systems show that tests wi 



16|, 



2l|. 



1 be satisfied by the small-Ci^2,3,4 



and the weak field PPN parameters conditions 

In IPN weak field experiments, the preferred frame PPN parameters are |a;i| < 10~^ 
and 1 02! < 10~^ jd, I44I. Putting both of them to zero leads to two conditions with two 



parameters 



C2 



{-2cl - C1C3 + c|)/(3ci), C4 = -c|/ci; 



• Ci3 = 0, Ci4 = 0. 

The latter case (ci4 = 0) will cause the propagation of the linearized wave modes of spin-0 
and spin-1 with infinite velocities ll|, [ij]. If it is shown that the difference between GR 



and ae-theory is extremely small b y fu ture 2PN light deflection experiments, it, together 

will improve the above conditions further 



with the IPN conditions given by 
(«! = a2 = plus Ci4 = 0): 



12 



C2 



-2ci/3, C3 — Ci, C4 — —Ci, 



• Ci3 = 0, Ci4 = 0, 

which implies that the velocities of spin-0 and spin-1 modes waves are infinite in both 
cases if there is no deviation between ae-theory and GR in weak field experiments. Other 
experiments focusing on 2PN periastron advances in the binary pulsars and the solar system 
and on preferred frame effects can provide more constraints on the action parameters, Ci^2,3,4- 
These experiments need 2PN equations of motion in ae-theory. The simplest case is the free 
fall of a test particle in the Sun's gravitational field. It associated metric in the isotropic 
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coordinate, which needs to extend goo in Eq. fll30p to is 



c^r 



2gMQ 2g^Ml /3 
^ " c*r^ ^ [ 2 



1 \ G'^Ml 

Cl4 



c'dt' 



+ 



2gM^ (3 1 



c^r 



12' c6r3 
(dV + rMn^ 



2 4 y cr 

where again only one non-GR parameter, Ci4, appears in goo and (7,^. 



(132) 



VII. CONCLUSIONS 



In this paper, we obtain the 2PN approximation of ae-thoery by Chandrasekhar's ap- 
proach. It shows a more comprehensive picture of the structure of 2PN approximation than 
scalar-tensor theories. Our works obtain five PPN parameters in IPN, and they are consis- 



tent with previous works 



I2,ll3|,ll4. Meanwhile, a 2PPN parameter in ae-theory, l, which 



is a non-GR parameter in the 2P N qp n, is discussed and compared with other theories, such 
as a multiscalar-tensor theory {25] 



30| and a scalar-tensor theory with an intermediate range 
force [35] . It is shown that l is zero in ae-theory, which means that ae-theory, similar to GR, 
does not have the third no linearity of Newtonian potential in goo. For future applications in 
deep space laser ranging missions, we derive the 2PN metric for light propagation in the case 
of nonspinning point masses as a simplified model for the solar system. The deviation of 
the resulting 2PN deflection angle of light between GR and ae-theory is dependent only on 
Ci4, which gives another constraint of action parameters, Ci^2,3,4- If future experiments show 
Ci4 is zero, it means that the linearized waves with the spin-0 and spin-1 modes in ae-theory 
will propagate with infinite velocities. 
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APPENDIX A: CONSTANTS AND 



There are some constants like Cy ^ in the coefficients of the metric, whose expressions 



are 



C 



(3) 



N 



,it 



c\ + 2ci — c{ 



, Ci3 / 1 ^2 1 

H -C13C23 + C3 - -C2C4 - C3C4 

C123 ^ 



(Al) 



1 - |C14 



7 9 3 



+ 7^ (C2 + C12 + C1234)' 

^Ci23 



/^'vij -1 



3 (Ci2 + C24)(C23 - C4) Ci4„i 



(A2) 
(A3) 



1 



4ci23(ci-c2 + 2ci) 



+ 6C2C3C14 + C3C4C234 + 4C4C3 + 3C1C4 



C1C4C23 - 2C1C4C12 - 4C1C4C123 - 2C3C23 



4ci23(ci - cf + 2ci) 



4C3C123 - 2CiCi4 - 24ciCi23 



+ C3C4 + I7C1C4C23 + 18CiCi234 - C1C4 



- 2C3C4 - C3C4C23 - 2C3C123 



(A4) 



(A5) 



(5) 



+ 2c\ + 14cf C2 + 8cf C3 - 14cf C4 



8c?23(ci - c? + 2ci^ 

+ 12c^C2 + I8C1C2C3 — 13clc2C4 — Scfcl — 23c^C3C4 + 5c^C4 

— 26C1C2C3 — 8C1C2C3C4 + 4C1C2C4 — 32C1C3 — 2C1C3C4 + 8C1C3C4 
— 12C2C3 — 30C2C3 + 5C2C3C4 + 4C2C3C4 — 18C3 + 7C3C4 + 3C3C4 
-32ci - 8OC1C2 - 96ciC3 + 12ciC4 - 48ciC2 - 124ciC2C3 - 4ciC2C4 

o 0000 

— 76C1C3 + 4C1C3C4 — 4C1C4 — 12C2C3 — 14C2C4 — 24C2C3 — 24C2C3C4 
-2C2C4 - 12cl - IOC3C4 - 2C3C4 ) , 



(A6) 



c 



(5) 

Ki 
NN.i 



8c?23(ci-cf + 2Ci) 



+ 184 + 28c^C2 + 46c^C3 + 32c^C4 + 10c?c^ 
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+58cfc2Cs + 51cfc2C4 + 50clcl + 73cfc3C4 — cfc^ + 2OC1C2C3 



,2 „2 



+20C1C2C4 + 52C1C2C3 + 80C1C2C3C4 + 34C1C3 + 48C1C3C4 + 10C2C3 



+2OC2C3C4 + 22C2C^ + 29C2C^C4 + I2C3 + 7C^C4 + C3C4 + 16c(c2 

+ 16c^C3 — 20c^C4 + I6C1C2 + 28C1C2C3 — I6C1C2C4 + I2C1C3 — I6C1C3C4 

-4C2C3 + 2C2C4 - 8C2C3 + 8C2C3C4 - 2C2C4 - 4c3 + 6C3C4 - 2C3C4 ) , (A7) 



(5) 



1 



+ 2CiCi34 + C?C234 + C3C234 



2ci23(ci-c2 + 2ci) 
+2C1C23C34 + 2C12C3C4 + 2CiC23 — 2CiC4 



(A8) 



^1 11 

2 - -Ci - C2 + -C3 - -C4 

2 - ci4 [ 4 4 8 

C12 + C24)(7ci + 2C2 + 3C3 + 5C4) 



4c 



123 



C4(Ci2 + C24)^ 



(A9) 



4-2ci 



5 1 

2C2 + ^C3 



15 
"8 



■C4 



(Cl2 + C24)(10Ci + 6C3 + 5C24) C4(Ci2 + C24)' 



4c 



123 



8c2 



(AlO) 



2 - Ci4 



c^ — cf + 2ci 



1 3 

2ci - -C24 - -C3 



(Ci + 2C3 - C4)(C3 + C1C4) - C4C13 ( -Ci3 



2r2 



1 



1 



/ 2 2 , o \ I o'^13C23 + C3 — -C2C4 — C3C4 

C123(C3 -cf + 2ci) V2 2 
2Ci3C23(ci3 + C3 - C4) 



C123(C3 - q + 2ci^ 



(All) 



(3) 



Ci4 

2 



Cl4 - C2 + 



(Ci2 + C24)(C23 - C4) 
2ci23 



c^ — c\-\- 2ci 



(Ci + 2C3 - C4) (C3 + C1C4) - C4C13 -Ci3 - 1 



^ c\^{c2Ca + 2C3C4 - C13C23 - 2C3) ^ 2ci3C23(ci3 + C3 - C4) 

Ci23(ci - + 2ci) 



Cl23(ci - Cf + 2ci) 



(A12) 
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In the Eq. (ITOl) . we introduce a variable A, which satisfies a Poisson's equation as 



^ (10C13 + 12c, - 3C4) + 



+ 



3 

2 »Ci23 
2C3 - 5C4 Ci2 + C24 / 1 , 1 



-123 



CIS 



4ci23 

C4 
2Ci23 
4Ci3 



iV.fcALfc + 



123 

Ci + 2C3 — C4 



2c 



123 



NtAN 



(3) 



(3) 

A^.fcAiTfe + 



Cl3 f 1 



LkAN^k - KkAN^k 

4Ci3(2-Ci4) 



7re;(Arat + aNt) - . ^ 
C123 C123 \2 



Cl3 



2c 



123 



,ttt 



7 



- + -ci4 NkN,kt + 

C123 V2 4 y 

Ci4 / 3 C1 2 + C24 
2Ci23 

C12 + C24 



C123 
N,ut 



2-r^^ + r^ + r^^+ 4C,23 



(2ci - 3c2 - C3 + 3C4) 



2-C14V2 4""^' 



+ 



Cl23(2 - C14) VS 



:Ci4 



. 7 9 3 ^ 

I + -C1 + -C2 + -C3 + 2C4 



(C2 + C12 + C1234, 



(A13) 
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